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Abstract 

A one-dimensional cellular automaton with a probabilistic evolution rule can generate 
stochastic surface growth in (1 + 1) dimensions. Two such discrete models of surface 
growth are constructed from a probabilistic cellular automaton which is known to show 
a transition from a active phase to a absorbing phase at a critical probability associ- 
ated with two particular components of the evolution rule. In one of these models, 
called model A in this paper, the surface growth is defined in terms of the evolving 
front of the cellular automaton on the space-time plane. In the other model, called 
model B, surface growth takes place by a solid-on-solid deposition process controlled 
by the cellular automaton configurations that appear in successive time-steps. Both 
the models show a depinning transition at the critical point of the generating cellular 
automaton. In addition, model B shows a kinetic roughening transition at this point. 
The characteristics of the surface width in these models are derived by scaling argu- 
ments from the critical properties of the generating cellular automaton and by Monte 
Carlo simulations. 
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1 Introduction 



Cellular automata are discrete mathematical systems that follow non-equilibrium dynam- 
ics [0. Typically a discrete variable with a finite number of possible states evolves in discrete 
time-steps, simultaneously at all points of a discrete space (for example, a lattice) on which 
it has been defined, following a local evolution rule. The simplest of these forms, known as 
elementary cellular automata, have two states per site and only nearest-neighbour interac- 
tions. The evolution of a <i-dimensional cellular automaton is by itself a process of surface 
growth in (d + l)-dimensions. Besides, it can generate other surface growth models, for ex- 
ample, by considering the configurations of a cellular automaton in successive time-steps as 
sieves with pores at the sites that have a particular state and through these pores particles 
are deposited on the corresponding sites of a substrate lattice. 

A cellular automaton with a probabilistic evolution rule will generate a stochastic growth 
process. When the noise due to the probability is of a non-thermal kind it may bring about 
a phase transition even in a one-dimensional system with short range interactions || . Such 
a phase transition in a cellular automaton will induce a corresponding effect in the surface 
growth models generated from it. Generally two kinds of phase transitions are observed in 
surface growth models : (1) morphological transitions such as kinetic roughening transitions, 
and (2) depinning transitions. While morphological transitions are due to a critical change 
in the rate of growth or deposition, depinning transitions occur by a critical increase in the 
force that drives the surface so that it overcomes the pinning force ||. Examples of such 
transitions in (1 + l)-dimensions are few, most of which are related to a Directed Percolation 
process : kinetic roughening transitions in (l + l)-dimensions have been observed in models of 
polynuclear growth H, growth by absorption and desorption of solid on solid []5| and fungal 
growth |Bj whereas depinning transitions in (1 + l)-dimensions were originally observed in 
models of quenched point disorder [|7j and porous media [§]. 

This paper reports a study of surface growth phenomena generated by a one-dimensional 
probabilistic cellular automaton that shows a phase transition in the universality class of 
Directed Percolation. 
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2 Two models with depinning transition 



The one-dimensional probabilistic cellular automaton (PCA) , which is the generator of two 
surface growth models studied in this paper, is defined as a line of sites with a binary variable 
dj G {0, 1} assigned to each site i. A site is said to be occupied if ai = 1 and unoccupied 
otherwise. The automaton evolves from any given initial configuration {a! " 1 } by updating 
the variables at all lattice sites simultaneously at each time-step. Each site follows a local 
rule of evolution that involves only nearest-neighbour interactions. The rule is specified by 
a set of eight components [afl^ af \ aj+J h- > af +1 ^ corresponding to the 2 3 distinct three-site 
neighbourhoods of interaction : 



t : 111 110 Oil 101 010 100 001 000 

t+ 1 : T T T T T ~o~ ^ 

■ ' 

1 with probability p 
with probability 1 — p 

The probabilistic behaviour enters the cellular automaton as a non-thermal noise added 
specifically to two mutually symmetric rule components. 

Based on the evolution of this PCA from disordered (uncorrelated) initial states two 
different models of surface growth in (1 + 1) dimensions are constructed : 

(1) Model A — The initial configuration of the PCA forms the substrate with the occupied 
lattice sites acting as growth centers. The propagating front of occupied sites on the space- 
time plane defines a surface. Here growth occurs in the direction of increasing time which 
therefore is the vertical direction along which the surface height will be measured. Overhangs 
may exist on the front but they are ignored so that the height profile is a single-valued 
function. The height of the surface hi(t) above a point i on the substrate at time t is given 
by the time coordinate of the highest occupied site in the i-th column which means that in 
this model all sites in a column in the space-time plane between the substrate site and the 
highest occupied site are assumed to be filled. This is a growth model of the propagating 
interface type. 

(2) Model B — This is a model of the deposition type. Surface growth takes place on a one- 
dimensional substrate (a line of empty sites) by depositing particles in successive time-steps 
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according to the configuration of the cellular automaton at those times. At time t a particle 
is added to the i-th column, i.e., the height of the z-th column is increased by unity if the 
corresponding site in the PCA is occupied at that time (of = 1). The height variable is 
defined by : 

fc(*) = E«S T) - ( 2 ) 

T=l 

A surface grown in this way does not have overhangs. Figure 1 shows typical examples of 
surfaces of models A and B generated by the evolving cellular automaton. To illustrate the 
procedure of growth in the two models outlined above the detailed mechanics are shown 
schematically for a small part of the lattice. 

The purpose of having disordered initial states of the PCA is to avoid the appearance of 
absorbing states at the beginning of its evolution. The PCA is known to have three absorbing 
states |l(J : the zero configuration (• • • 0000 • • •) and two mutually symmetric configurations 



given by sequences of alternating and 1 (• ■ -0101 • • -, • ■ ■ 1010 ■ • •). Occurence of these 
absorbing states leads to either no growth (for the zero configuration) or a trivial growth 
process (for the other two absorbing states). A disordered state is constructed by assigning 
each lattice site the value 1 with a predefined probability and otherwise, independently of 
the values at other lattice sites; when this probability is | the state is said to be completely 
random. In general, any initial state other than the three absorbing states will serve this 
purpose. 

The growth models are characterised by the properties of their surface width which, for 
a substrate of linear size L, is defined in the usual way || : 



w 2 (L,t) = ( [hi(t)-h(t)] 2 j , (3) 

h(t) = j- J2i=i hi(t) is the average height of the surface above the substrate at time t. The bar 
denotes average over the system and the angular brackets (• • ■) denote configurational 
average. In the addition to the width, a moving surface is characterised by its average 
velocity of propagation, defined by || : 

vt = §jKt) ■ (4) 

For the PCA model there exists a critical value of p, observed to be p c ~ 0.75 0, at which 
a continuous phase transition occurs with respect to the density of occupied sites in the steady 
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state; Monte Carlo simulations JTIJ of the PCA indicate that the phase transition belongs 
to the universality class of Directed Percolation. In the absorbing phase (p < p c ) all initial 
configurations other than the two mutually symmetric absorbing states following the rule 
in equation ([3]) evolve to the zero configuration (all sites unoccupied) which is a absorbing 
state of the PCA. Evolution of the cellular automaton stops when the zero configuration 
is reached. The presence of occupied sites in the cellular automaton configuration at any 
time-step produces a moving surface in the growth models. The appearance of the zero 
configuration results in pinning of the surface []. In the active phase (p > p c ) all initial 
configurations other than the three absorbing states evolve to a active steady state with a 
non-zero density of occupied sites. This results in unpinned surfaces in the growth models. 

Thus the critical point p c of the cellular automaton marks a depinning transition in the 
growth models, from a pinned phase for p < p c to a moving phase for p > p c . The order 
parameter for the phase transition in the cellular automaton is the density of occupied sites 
Poo in the steady state which in the supercritical region decreases continuously to zero as the 
critical point is approached || : 

Poo OC (p - p c Y, P~*Pc ■ (5) 

The appropriate order parameter in the growth models in general is the asymptotic velocity 
foo of the surface in the moving phase || : 

Woo OC (p - Pcf, P-^Pc ■ (6) 

The following sections report the investigations on the depinning transition in models A 
and B by scaling arguments and Monte-Carlo simulations. 

3 Properties of Model A 

By definition Model A is the direct evolution of the PCA from disordered initial states 

viewed as a growth process that excludes surface overhangs. Therefore the properties of the 

surface near the depinning transition can be derived from those of the critical properties 

of the PCA. Figure 2(a,b) shows typical examples of evolution of the surface width w(L,t) 

1 Though pinning centers and forces are not explicit in these models, pinning of a surface growing in 
(d + l)-dimensions can be thought of as a d-dimensional system entering an absorbing state where there are 
no particles. 
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with time for substrates of different sizes L obtained by Monte Carlo simulation of Model 
A. As it is usual for correlated growth on substrates of finite size the surface width gets 
saturated after a crossover time t x . Both above and below the depinning transition point 
p c the saturated value of the width, denoted by w(L, oo), shows two regimes separated by 
a crossover length L x — (1) for L <C L x it increases with L indicating the appearance of a 
'rough' (or, L-dependent) regime; (2) for L>L X , w(L, oo) tends to saturate to a constant 
value as L — > oo which indicates a 'smooth' (or, L-independent) regime. The values of L x 
and w(oo, oo) depend only on the value of the noise parameter p. Therefore the saturated 
width w(L, oo) is expected to follow a scaling relation of the form 

w(L,oo) ( L \ 

—7 r oc /<> — , (7) 

1(7(00,00) \-^x/ 

where /<, /> denotes the scaling function in the pinned (p < p c ) and moving (p > p c ) phases 
respectively. 

The crossover length L x provides a characteristic lengthscale of the model. Because of 
the scaling hypothesis that the model has a unique diverging lengthscale at the critical point, 
L x must diverge as the correlation length of the PCA as p approaches p c from above and 
below : 

L x oc \p-p c \- v± , (8) 

where Vj_ is the critical exponent for the correlation length of the PCA. Since the growth 
of the surface occurs in the direction of time, the saturated surface width is a measure of 
temporal RMS fluctuations in the steady state. Consequently, in the limit of an infinite 
substrate, the quantity w(oo, oo) will diverge as the correlation time of the PCA : 

w(oo,oo) oc \p-p c \~ n , (9) 
This reduces the scaling law (|7|) to the standard form of finite size scaling [JTTJ] : 



w(L, oo) \p - p c p oc /<>(£ \p - p c r ) . (10) 



Since the PCA is believed to be in the universality class of Directed Percolation [jnj, the 
critical exponents for L x and w(oo, oo) are expected to be v± = z/^ p w 1.097 and v\\ = 
uP p w 1.734 respectively O. 
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At p c the crossover length L x diverges and only the rough regime exists. The surface 
width is then expected to follow the Family- Vicsek scaling law [[0| : 



w(L,t)xL a f(J^j, p = Pc . (11) 

The growth exponent A and the roughness exponent a are defined as w(L,t) oc t x , t <C t x 
and w(L,t) oc L a , t^>t x respectively and the dynamic exponent z occuring in ( |TTD is given 
by t x oc L z , z — a/A. The value of the exponents a, A and z at p c can be obtained by using 
the scaling hypothesis : corresponding to a uniquely diverging lengthscale of the model there 
exists a uniquely diverging timescale. The saturated width w (L, oo) in this model measures 
temporal RMS fluctuations in the steady state and in the thermodynamic limit {L — > oo) it 
must diverge as the correlation time as p approaches p c . For a finite-sized system (substrate) 
at p = pc, w(L, oo) oc L v \\l vx which gives the roughness exponent : 

<W = — • (12) 

Similarly the crossover time t x is another way of defining the characteristic timescale for the 
model and shall have the form t x oc L u ^ u± at p = p c . Therefore the dynamic exponent is 
given by : 

z = ^. (13) 
v± 

Since the phase transition of PCA is believed to be in the universality class of Directed 



Percolation |10] 



a p=Pc = z = -jjp w 1.58 . (14) 
Therefore the growth exponent at p = p c becomes equal to unity : 

X p=Pc = 33* = 1 . (15) 

z 

To confirm the results obtained above from scaling considerations Monte Carlo simula- 
tions of the growth process in Model A were performed on substrates of sizes L = 100 to 
L = 10000. A disordered substrate is constructed by assigning to each lattice site the value 
1 with probability | and otherwise. For each substrate size 2000 independent simulations 
were done for the purpose of averaging out statistical fluctuations. 
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Figure 3 shows the results of applying scaling relation ( |TD| ) to the simulation data for the 
surface width. Data for w(L, oo) vs L in the pinned phase, for different values of p close to 
p c , collapse to a single curve for the following values of the transition point and the critical 
exponents : 

p c = 0.7513 ±0.0004; (16) 



u± = 1.1 ±0.01, v\\ = 1.73 ± 0.01 . (17) 

Similarly for the moving phase, values of the critical point and exponents were obtained by 
collapsing the simulation data for w(L, oo) vs L for different values of p close to p c : 

p c = 0.7511 ±0.0004; (18) 



u± = 1.1 ±0.01, v\\ = 1.74 ±0.01, (19) 

These are found to agree, within the limits of error, with the corresponding values for the 
pinned phase. Error bars indicated here and elsewhere in the paper are statistical errors. 

The estimate of the depinning transition point p c obtained here agrees well with the value 
of the critical point of the PCA measured by using defect dynamics ||10|| . The estimates of 
the critical exponents are in good agreement with the values for the universality class of 
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To measure the surface exponents for p = p c simulations were performed at p = 0.7512. 
The time taken for the PCA to reach a steady state is very large and the width of the surface 
could not be measured till saturation. Therefore a direct determination of a and z at p c 
was not possible with the available computer facility. However the surface width (averaged 
over 5000 intial configurations) was observed to grow linearly on a substrate of size L = 10 4 
for 10 5 times-steps (Figure 4), thus establishing X p=Pc = 1 and a p=Pc = z, in agreement 
with equation fllBl) . Analyses of simulation data far from p c show that the 'rough' regime 
occurs for only very small system sizes (L <C L x ) and the dependence of w(L, oo) on L 
is logarithmic for the most part of the crossover region. The 'rough' regime grows as p 
approaches p c and w(L, oo) becomes entirely a power-function of L at p = p c . Thus a true 
'rough' surface (in the sense of a > 0) exists only at the transition point. 
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The last result of this section is concerned with the velocity of the surface in the moving 
phase. Simulations of the growth process A for different values of p show that the surface 
in the moving phase attains the maximum velocity asymptotically, = 1, irrespective 
of the value of p (Figure 5). Consequently cannot serve as the order parameter for the 
depinning transition in Model A. The appropriate order parameter is the asymptotic density 
of columns with the maximum height which is equal to the density of occupied sites in the 
steady state of the generating PCA and hence it decays by the power-law @ with the critical 
exponent (3 = f3 DP « 0.276. 

4 Properties of Model B 

Model B defines a deposition process of the solid-on-solid type, i.e., the bulk of the aggregate 
is compact (no vacancies) and there are no surface overhangs. Consider first the moving phase 
(p > p c ). Since there is no interaction between the columns of the aggregate, correlations do 
not develop along its surface; hence the surface width in the moving phase does not saturate. 
When the PCA evolves from a disordered initial state the process can be described by a 
continuum equation similar to the one describing a ordinary random deposition process || : 

^h(x,t) = F + ri(x,t), (20) 

where F is the average number of particles arriving at a site and rj is the noise term (with 
zero configurational average : (r](x,t)) = 0) that describes the fluctuations in the deposition 
process. The difference between the random deposition model and this model lies in the noise 
correlations. While the noise in random deposition is uncorrelated, the spatial and temporal 
correlations developing in the generating PCA appear as noise correlations in model B. Since 
the initial state of the PCA is disordered the growth process in the first time-step looks like 
random deposition. But correlations in space and time develope in the PCA configuration 
as it evolves owing to the nearest-neighbour interactions defined in the evolution rule ([J). 
For values of p away from p c , the correlation length £ and the correlation time r of the PCA 
are finite which means the noise in the deposition process is correlated over a short range 
described by exponential decay of the correlation function [[| : 

(ti(x, t)r](x\ £')> ~ e-\ x - x '\'t e'^'^. (21) 
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The effects of the correlated noise is prominent for times less than r. Since deposition at 
time-step t + 1 is allowed only at or next to sites where deposition has taken place in the 
previous time-step t [along with the restrictions set by the PCA evolution rule ([!])] the 
height fluctuations about the mean value are larger and hence the growth of the surface 
width is faster than that for random deposition. For times greater than r the noise appears 
uncorrelated and the growth exponent A will reduce to the random deposition value X RD = |. 
As p approaches p c , £ and r increase and it takes longer for the growth process to reach the 
random deposition limit. This is illustrated in Figure 6. Finally, at p = p c , both £ and r 
diverge and the correlations in the noise are long-ranged, represented by a power-law decay 
of the noise correlation function |5| : 

{r]{x,t)r]{x' ,t')) ~ \x - x'\~ 2P/v± I* " ■ ( 22 ) 

where /3, v± and vn are the critical exponents for the order parameter, correlation length and 
correlation time of the PCA respectively. The asymptotic value of the growth exponent A 
at p c can be derived from the continuum equation fl20|) where rj(x,t) follows the correlation 
function (|22| ) instead of the delta-correlated noise of a random deposition process. The 
surface width for a infinite substrate is found to increase as the following power function of 
the time : 



w 



(oo,t) oct 1 "^". (23) 



Using the indication that the phase transition in the PCA belongs to the universality class 
of Directed Percolation |nj, the value of the growth exponent at p c is expected to be : 

X P=Pc = 1 - ^pp « 0.841. (24) 

At this point the surface grown on a finite substrate enters the pinned phase. 

The behaviour of the velocity of the surface in the moving phase can be analysed as 
follows. The average height of the surface above the substrate at time-step t is equal to the 
number of particles deposited per site in the interval [0, t] averaged over the entire system. 
The velocity of propagation v t of the surface at time-step t, defined by Eq. (||), is therefore 
equal to the average number of particles deposited per site in the t—th time-step. This 
number is precisely the density of occupied sites p t of the generating PCA at time t : 
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v t = Pt- (25) 

In the steady state, as p appoaches p c , the asymptotic velocity of the surface will thus 
decrease to zero by the power- law followed by p^. This implies, after comparing equations 
(|5|) and ©, that 9 = (3. Since the phase transition of the PCA is known to be in the 
universality class of Directed Percolation |H| , the critical velocity exponent is given by : 



= [3 DP « 0.276. (26) 

The model B in the moving phase was studied numerically by Monte Carlo simulations 
of the growth process on a one-dimensional lattice of L = 10 4 sites for various values of p 
(p > p c ) and results for each value of p were averaged over 10 3 independent realisations of 
the process obtained by constructing the disordered initial states of the PCA independently 
of one another. Figure 6(a) shows the surface width as a function of time for different values 
of p, drawn on double-logarithmic scale. To determine the growth exponent A an effective 
exponent X t is defined as the local slope of the curves drawn in Figure 6(a) : 

^ = log[w(L,t)/w(L,t/b)] ^ 
log b 

With 6 = 5, Xt was calculated for different values of p upto t = 10 5 time-steps and plotted 
against 1/t (Figure 6b). All curves for p > p c tend to the random deposition value X RD = |; 
though the curves for p close to p c do not reach this value in 10 5 time-steps they clearly show 
a tendency to do so at a larger t. For p = 0.7512(^ p c ), in order to remove finite-size effects, 
the growth exponent was obtained by extrapolating the maximum values of the effective 
exponent measured for different substrate sizes : 

X p=Pc w max [X t {p = 0.7512)] = 0.837 ± 0.011, (28) 

which is close to the expected value for an infinite substrate [equation (0)]. In the curve for 
p = 0.7512 in Figure 6(b) a decrease observed in the value of X t after it reaches a maximum 
is due to the finite size of the substrate used in the simulation and it shows the beginning 
of a crossover to saturation of the surface width, the system having just entered the pinned 
phase at p c . 

The asymptotic velocity of the surface obtained from the simulation data (Figure 7) 
show that the power-law behaviour of v oo as p — > p+ [equation (H)] is best satisfied with : 
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p c = 0.7514 ±0.0004. 



(29) 



This independent estimate of p c agrees well previous estimates ]TD| and those obtained from 
the study of model A in this paper. The value of the velocity exponent was measured to be 



9 = 0.274 ±0.007, (30) 

which agrees, within the limits of error, with the expected value given in equation (p6|). 

Consider now the pinned phase. For p < p c growth of the aggregate stops in finite time 
(which depends on the parameter p) due to the appearence of the absorbing state of the 
generating PCA. The time for which growth takes place is equal to the transient time T of 
the generating PCA, i.e., the number of time-steps required by the PCA to evolve from a 
disordered initial configuration to the absorbing zero configuration. Like that of model A, 
the surface width w(L, t) saturates after a crossover time t x for a substrate of finite size. The 
behaviour of the saturated width w(L, oo) too resembles that of model A : after an initial 
Independent regime it passes over to a L-independent regime (Figure 2c) beyond a crossover 
length L x . Thus for p < p c the surface in model B is smooth (in the sense of a = 0) in 
the thermodynamic limit (L — > oo). Only at the depinning transition point p c , where L x 
diverges, the surface is rough (a > 0) for all substrate sizes and the Family- Vicsek law ([11]) 
is expected to hold true. Since the saturation of the surface width is forced by the PCA, the 
crossover time t x will be the same as that of model A which means the dynamic exponent 
at p c is given by the Directed Percolation value : 

t x (p = p c ) oc L z , z = z DP « 1.58. (31) 

At p c , just before the surface enters the moving phase, the roughness exponent is given 
by : 

a p=Pc = z ■ \ p=Pc = z DP {\ - (3 DP /vf P ) « 1.33. (32) 

This result for a could not be verified by Monte Carlo simulations on available computer 
facility as it takes enormous times for the surface width on substrates of reasonable sizes to 
reach saturation at p c - 
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Since the surface width does not saturate in the moving phase, it becomes infinitely 
large with time. Therefore the depinning transition point also marks a kinetic roughening 
transition from a smooth surface in the pinned phase to a infinitely rough surface in the 
moving phase. 



5 Discussion 

In the two models studied in this paper the dynamics of surface growth are derived from 
the evolution rule of a one- dimensional probabilistic cellular automaton. While model A 
presents the growth of an aggregate of interacting particles, model B presents the growth 
of an aggregate of non-interacting particles deposited by a correlated mechanism. This 
results in the formation of a correlated surface in model A and an uncorrelated surface in 
model B. The correlated mechanics of deposition in model B leads to growth of the surface 
width of the aggregate at the depinning transition point p c with an exponent A larger than 
that offered the ordinary random deposition model. This result follows the trend that the 
growth exponent increases in the presence of correlations in the deposition process, observed 



previously in ballistic deposition models Also the surface properties of the models at 
the depinning transition are all related to the critical properties of the generating PCA. It 
is yet to be seen the effect of introducing in these models the additional dynamics of surface 
restructuring. 

Finally, a comment on the scaling relation : 



9 = (z- a)v, (33) 

which relates the exponents for roughness, growth and velocity of the surface to the exponent 
v for the correlation length along the surface at a depinning transition. For model A at p c , 
v = z/j_ and a = z so that 6 = 0; this is true as the velocity in the moving phase is 
independent of the level of noise p which provides the driving force. For model B at p c , 
the values of 9, z and a from equations (p6|) , (^TJ) and (0) give v = v\_\ this is not true. 
Since there is no inter-columnar interaction at the surface of the aggregate, correlations do 
not develop along the surface and an exponent v cannot be associated with it; the scaling 
relation (^) therefore does not hold in this case. 
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Figure Captions 



Figure 1. (a) An example of the evolution of the probabilistic cellular automaton following 
the rule in equation (1) for p = 0.8. Sites (i, t) on the space-time plane with automaton 
value a® = 1 are marked dark while those with automaton value af^ = are left blank. 

(b) The surface of model A and (c) that of model B at the 250-th time-step of the 
automaton evolution shown in (a), (d) Schematic diagram showing the construction 
of the surfaces of models A and B : model A is constructed from the PCA pattern by 
filling all vacant sites below the highest occupied site in every column whereas model 
B is constructed by considering the occupied sites in the PCA pattern as particles and 
dropping them one above the other in a ordered way in every column. The surfaces 
are shown by lines joining the highest point in each column. 

Figure 2. The development of the surface width in time on substrates of different sizes : 
(a) model A in the pinned phase (p = 0.7), (b) model A in the moving phase (p = 0.8), 

(c) model B in the pinned phase (p = 0.73). 

Figure 3. Data collapse obtained by applying the finite-size scaling relation in equation (10) 
to the simulation data for the saturated surface width of model A in the (a) pinned 
and (b) moving phases. 

Figure 4. The growth of the surface width of model A at the depinning transition point 
p c ^ 0.7512. The solid line shows the result of simulation on a lattice of 10 4 sites; 
the dashed line represents a linear growth of the surface width : w(L,t) oc t. At large 
times the surface width is seen to grow linearly with time as it becomes parallel to the 
dashed line. 

Figure 5. The average velocity of the surface, defined by equation (4), in the moving phase 
of model A is plotted against time. For all values of p, a few of which are shown here, 
the velocity asymptotically attains the maximum possible value of one lattice constant 
per time-step. 

Figure 6. (a) The growth of the surface width of model B in time for different values of p 
in the moving phase. Simulation data are shown for the following values of p : 0.7512 
(top), 0.754, 0.756, 0.758, 0.76, 0.77, 0.78, 0.79, 0.8 and 0.9 (bottom). The curves do 
not appear to saturate even after 10 5 time-steps of evolution, which agrees with the 
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theory, (b) The effective growth exponents At derived as local slopes of the curves 
shown in (a). The curves for p = 0.754, . . . , 0.9 tend to the random deposition value 
A = 0.5ast-^oo. The downward bending of the curve for p = 0.7512 is a signature of 
relaxation toward a saturated state; this is not easily evident from the corresponding 
curve in (a). 

Figure 7. The asymptotic velocity of model B in the moving phase for different values p 
near p c . The inset shows that it decreases to zero by a power-law as p approaches p c . 
The slope of the dashed line gives the value of the velocity exponent 9 (see text). 
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Figure 2(c) 



10 - 



1 - 



L = 100 
L = 200 
L = 500 
L= 1000 
L = 2000 
L = 5000 








+ + 



2000 4000 6000 8000 10000 

L 



0.1 



10 



100 



1000 



10000 



1.2 



Figure 3(a) 



1 - 



p = 0.735 
p = 0.74 
p = 0.742 



+ 

□ 

m 



□ 



• + 



0.8 - 



»□ 



□ 



HQ 



+ 



0.6 h + + 

/ 

ft 

U □ 

3 



0.4 



8 
_T 

5 



+ 



0.2 h 

□ 

I 



3000 
2500 
2000 
1500 
1000 
500 




□ 



□ 



i □ 5 5 if 



□ 



□ 



□ 
□ 

J. 



+ 



+ + 



+ + + + + 



f 



2000 4000 6000 

L 



8000 10000 











20 



40 60 

L |p - 0.7513| 



80 



100 



120 



1.1 



0.1 



Figure 3(b) 
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